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I. INTRODUCTION 



The purpose of this Letter is to discuss the relative motion of a spinning test particle 
in a gravitational field with respect to a nearby geodesic observer. To compare the theory 
with observation, it proves useful to express such relative motion in a quasi-inertial Fermi 
normal coordinate system that can be set up along the worldline of the reference observer 
[lj]. This approach has been recently employed for the general motion of a test particle; in 
particular, astrophysically significant results have been obtained for ultrarelativistic relative 
motion, i.e. motion with relative speed above a critical speed given by c/y/2 ~ 0.7c p. The 
present work aims to extend previous results to the case of a spinning test mass. We choose 
units such that c = 1 in the rest of this Letter. 

The motion of an extended body in a gravitational field can be described using the 
Mathisson-Papapetrou-Dixon equations jjj. If we neglect the quadrupole and higher mo- 
ments of the particle, then the Mathisson-Papapetrou ( "pole-dipole" ) equations suffice; these 
are given by 

npc 1 

dr ~ 2 ^ ' [ } 

DS^ 

—— = p» u v - P v u». (2) 
dr 

Here r is the proper time, = dx^/dr is the four- velocity, is the four-momentum of 
the particle and S^ u is its spin tensor. Suppose that the particle has mass m and spin so 
and is moving in the gravitational field of an astronomical mass M 3> m. In general, the 
dipole interaction must be much smaller than the monopole interaction in our approximation 
scheme; thus, the spin-curvature force in equation ~ (GM/r 3 )so, must be much smaller 
than the Newtonian force (GMm/r 2 ). This means that the M0ller radius 4| of the particle, 
p := Sq/ttl, should be much smaller than the distance r between m and M. 

Equations ((H) and (J2J must be supplemented with additional constraints on the spin 
tensor. For extended bodies, the appropriate condition turns out to be 
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S^P V = 0. (3) 

In section HH we develop an iterative scheme for the solution of equations Section 
TTTlis devoted to a discussion of the spin-curvature force. The equations of relative motion 
are developed in section |IV| and specialized to the exterior Kerr spacetime in section El 



The effects of spin-curvature coupling along the rotation axis of a Kerr black hole are of 
particular relevance to current theoretical speculation on the nature of gamma ray bursts. 
We conclude with a brief discussion of our results in section I VII 



II. APPROXIMATION SCHEME 

It is necessary to develop an iterative approximation method for the solution of equations 
(01-0. We note that these equations imply that u^DP^/dr = and 

-SfMuS^ = si, (4) 

where So is the constant magnitude of the spin of the particle. Differentiating (jSJ) with 
respect to proper time and using (j2J) results in 

DP 

(P ■ m)P^ - PV + S^—r^ = 0. (5) 

dr 

Next, we multiply (JHJ) by DP^/dr and use the antisymmetry of the spin tensor to get 

DP 2 

(P.u)_ = 0. (6) 

Thus either P • u = or P 2 = — m 2 , where m has the interpretation of the constant mass of 
the spinning particle. The signature of the spacetime metric is assumed to be +2 throughout 
this work. The multipole approximation method discussed in the previous section implies 
that P^ should be approximately parallel to w M except for small corrections due to the spin 
of the particle; in fact, multiplication of (J2J) with u v results in 

P" = -(p. u )u»-u v ^— . (7) 

dr 

We therefore assume that P • u ^ and write 

pv = mu^ + E 1 *, (8) 

where the extra term P M is the small contribution of the spin to the canonical momentum. 
From (|Sji . we find m = —P ■ u + E -u; on the other hand, multiplying (JHJ) with results in 

DP 

m 2 = (P-u) 2 + u^^. (9) 

Hi I 
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It thus follows from equations © and (jHJ) that E ■ u = mil — y/1 — e), where 

1 DP 

e = —u^—^ (10) 
m z dr 

is ~ (GM/r)(p/r) 2 < 1. 

Differentiating (jHJ) with respect to r and using (JTJ) we get 

m™-- l -R» u»S°*-^- fill 



To determine E 11 , we substitute (jHJ) in (0) to get 



1 „,.„DP„ 



yJT^l E»=(E- u)u» + — S^— A (12) 

m ctr 

It is clear from equations (jlU)) - (jl2|) that E^ is generally small and of order me; moreover, 
DE^/dr is a force that is of second order in the spin of the particle. In this way, it is 
possible to develop an iterative scheme based on the small parameter p/r <C 1. Restricting 
our treatment to effects that are of first order in the spin, we drop DE^/dr in (|TT|) and 
focus on the Mathisson-Papapetrou spin-curvature force 

r = -\r^s^. (is) 

It is useful to define the spin vector S 1 * of the extended test particle in general as 



S, = 7^V»^P"S p °, (14) 

where r/a^s = y/^g is the Levi-Civita tensor and e a ^ is the alternating symbol with 
eoi23 = 1- Equation (fTlj) implies that S^S^ = Sgj moreover, 

mS aP = r) a ^ 5 P 7 S s . (15) 

It follows from the results of this section that to linear order in spin, one can replace P M by 
mw M in equations (pPfj) and (fTK|) . 

III. SPIN-CURVATURE FORCE 

Consider the linear post-Newtonian gravitational field of a source with mass M and 
angular momentum J given by 

- ds 2 = -(1 + 2$)dt 2 - 4(A ■ dx)dt + (1 - 2$)5 ij dx i dx j , (16) 



where $ = —GM/r is the Newtonian gravitoelectric potential, A = GJ x x/r 3 is the 
gravitomagnetic vector potential and r = |x|. The gravitomagnetic field is given by B = 
V x A in analogy with electrodynamics. If an ideal torque-free test gyroscope is placed at 
x, its precession frequency is given by 

B = ^[3(x-J)x-Jr 2 ], (17) 



so that the gravitomagnetic field has the interpretation of a precession frequency in con- 
formity with the gravitational Larmor theorem. Expressing the metric tensor in (fTKJ) as 
9pu = Vnv + h^ u , we note that h 0Q = —2$, hij = —2§5ij and h 0i = —2Ai. The corresponding 
Riemann curvature tensor is given by 

Rpupa ~^{hpa,vp hup, pa hva,pp ^-/jp^cr)' (1^) 

Imagine a free test particle S initially at rest in the stationary gravitational field of the 
source. The particle carries an orthonormal tetrad frame A'Vn, where A M , Q s = is its 
local temporal axis and i = 1,2,3, are its local spatial axes. The curvature tensor as 
measured by S is given by the projection of (jl8)l on its local tetrad frame. To calculate this 
curvature in the linear approximation in the gravitational potentials, one may set = 5^ a \ 
then, the measured components of the Riemann tensor can be expressed as a 6 x 6 matrix 



£ B 
B -E 



(19) 



with indices that range over {01, 02, 03, 23, 31, 12}. Here £ and B are the electric and mag- 
netic components of the curvature, respectively, and are 3x3 symmetric and traceless 
matrices given by 

Eii = ^(<%-3x^), (20) 
C J 

Bij = -3— [x\P +x j J i + (5 il -5x i x :i )x- J], (21) 

where x = x/r and J = J/J. 

The spin-curvature force experienced by the extended test particle S is given by = 
f^\p a \ so that using equation |T3|) . 

/(°) = o, /W = B ij Sj , (22) 



where S( a ) = (0,s). More explicitly, this force can be written as 

f = -V(s-B), (23) 

which is the gravitational analog of the Stern-Gerlach force. 

As an illustration, let us suppose that S is initially at rest on the z axis with its spin 
along the z direction, which we take to be the rotation axis of the source; then, 

6G Js „ 



z 4 



(24) 



where z 3> 2GM. This Mathisson-Papapetrou force is repulsive (attractive) if S spins in 
the same (opposite) sense as the rotation of the central source. Moreover, it is interesting 
to note that this force is always smaller than the Newtonian force of attraction (GMm/z 2 ). 
Indeed, in realistic astrophysical situations, the Mathisson-Papapetrou dipole force would 
be very small in comparison with the Newtonian monopole force. This holds in general 
relativistic situations as well, as demonstrated in section 

If instead of being at rest, S is initially boosted with speed vq along the z axis, then the 
spin-curvature force experienced by the boosted particle can be obtained via a Lorentz 
transformation. In fact, the spin tensor remains invariant under a boost along the z axis, 
so that S 12 = —S 21 = So are the only nonzero components of the spin tensor. As discussed 
in the next section, the curvature tensor remains invariant as well. In any case, the nonzero 
components of F^ a > can be simply calculated to be 

F^= lo v f, F® = 7o f, (25) 



where 7 = (1 — Vq)~ 1 / 2 and / = 6GJs / ' z 4 is given by equation 

It should be emphasized that equations (J22j) - (|2l)j) are valid only to first order in the 
spin of the extended test particle. Similar results have been derived before using the other 
principal interpretation of the Mathisson-Papapetrou equations, namely, within the context 
of a classical point particle with "intrinsic" spin; that is, a point gyro that satisfies instead 
of © the Pirani supplementary conditions S^ v u v = |. 

IV. RELATIVE MOTION IN FERMI COORDINATES 

Imagine a free test "observer" O following a timelike geodesic in the exterior gravitational 
field of an astronomical source. Let A M ^ be an orthonormal tetrad frame that is parallel 



transported along the worldline of O. A Fermi normal coordinate system (T, X) can be set 
up in the neighborhood of this worldline based on A'Vv as the local axes such that O remains 
at the spatial origin of this coordinate system. Thus O has Fermi coordinates (T, 0) , where 
T is the proper time of O. We need to express the motion of a spinning test particle S 

™ = #, A ^-±.R^S^ (26) 

with respect to O. The equation of relative motion has been derived in general in []| and is 
given by 

^ + (o-r° e n^^'-^n (27) 

where the four-velocity w M in Fermi coordinates is r(l, V) with V = dT/dr. The requirement 
that the worldline of S be timelike can be expressed as 

r~ 2 = -goo ~ 2g 0i V l - g tJ V l V> > 0. (28) 

The metric in Fermi coordinates is given by 

g 0i = ~ F Ro Tlk (T)X^X k + ---, 

9H = Sij - i F R lk]l (T)X k X l + ■■■, (29) 



F Ra/3-ts(T) - R pvpa \^ ,.\ V , .X 7 r s y (30) 



where 

' ' ' \vpa* ( a )^ (/3) A ( 7 ) A {&)■ 

The Fermi coordinates are admissible in a cylindrical region with |X| < 1Z along the reference 
worldline such that 7Z(T) is a certain minimum radius of curvature of spacetime. 

In the next section, equation ()27|) is studied for the case of motion along the rotation axis 
of a Kerr source. The motion of spinning test particles in the Kerr field has been the subject 




of previous investigations jsj, |9j, llO, [ill, ll2j, ll3j] ; however, in the present work we consider the 
deviation of such motion relative to a reference geodesic in the context of a Fermi normal 
coordinate system that is constructed along the reference geodesic. This approach turns 
out to be of direct observational relevance for the motion of clumps in astrophysical jets as 
discussed in section EH 
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V. MOTION IN THE KERR FIELD 



The purpose of this section is to study the motion of S relative to O along the rotation 
axis of a Kerr source. The Kerr metric is given by 

- ds 2 = -dt 2 + £ (J^-dr 2 + dO 2 ^ + (r 2 + a 2 ) sin 2 9d0 2 

+2GM^(dt-asm 2 6d(j)) 2 (31) 

in Boyer-Lindquist coordinates, where S = r 2 + a 2 cos 2 9 and A = r 2 — 2GMr + a 2 . Here 
M and a > are respectively the mass and the specific angular momentum (J/M) of the 
source. We assume that the free reference particle O moves along the rotation axis on an 
escape trajectory. The geodesic equations of motion of O reduce to 



dt _ r 2 + a 2 dr / ~ 2GMr 

~ds~ 1 r 2 -2GMr + a? 1 ~d~s ~ V 7 + r 2 + a 2 ' ^ ^ 

Here 7 > 1 is a constant of integration such that for r — > 00,7 is the Lorentz factor of 

the particle as measured by the static inertial observers at spatial infinity. We integrate 

system with the initial conditions that at s = 0, t = and r = r > \^3a. Moreover, 

the spinning test particle S also starts from this same event with speed Vq > relative 

to O and moves along the rotation axis. To describe this relative motion, we establish a 

Fermi coordinate system along the worldline of O. In (t, r, 6, <fi) coordinates, the orthonormal 

tetrad frame % is such that 

(a) 

A% = (*, r, 0, 0), X\ 3) = {y-Hf, 7, 0, 0), (33) 

where t = dt/ds and r = dr/ds are given by (|32jl . The axial symmetry about the rotation 
axis implies that there is a simple rotational degeneracy in the choice of and ^( 2 )- The 
projection of the curvature tensor on the tetrad frame of O, given by equation (J3U)) . turns 
out to be independent of the explicit choice for and A m , 2 n and may be expressed as 
before in terms of £ and B such that £/k — diag (— |, — |, l) and B/q = diag (— |, — |, l), 
where 

r r(r 2 - 3a 2 ) - o Of a 3r ' - ^ 
(r 2 + a 2 ) 3 ' (r 2 + a 2 

Thus for a given r, the curvature measured by O is completely independent of 7. This 
remarkable fact is a consequence of the degenerate nature of the Kerr solution, which im- 
plies that the rotation axis of the exterior Kerr spacetime corresponds to two special tidal 



k = -2GM l„> , q = 2GMaj-^—^. (34) 
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FIG. 1: Plot of f versus T j {GM) based on the integration of equation for ro = 5 GM, p = 2, 
a = GM and T = 100, 200, 400, 600, 800 and 1000. 

directions for ingoing and outgoing trajectories; in fact, £ and B are invariant under a boost 
along the rotation axis 

The equation of motion of S relative to O is given by equations and ()27j) . where we 
need to evaluate the force term involving i? M VQ ,^(T, X) and S a ^(T, X) in the Fermi system. 
To simplify matters, we assume that R^^iT, X) « i? M ya/3 (T, 0), which can be determined 
by equation (}3"4")) . and S al3 (T, X) « S al3 (T, 0). As in the previous section, the spin of the test 
mass S is assumed to be along the axis of rotation, so that the only nonzero components 
of the spin tensor are S 12 = —S 21 = Sq. It follows that ^4° = qpTV and A 3 = qpT are the 
only nonzero components of A^, where T" 1 = yl — Z 2 + kZ 2 and p = So/m is the M0ller 
radius of S. With these simplifications and taking only the terms given in (|29~j) into account, 
equation (|27j) reduces to 



Z + k(l- 2Z 2 )Z = qp(l - Z 2 )T-\ 



(35) 



This equation takes a dimensionless form if all lengths are expressed in units of GM and p 
is replaced by p = p/(GM). Then, [GM/r)p = so/(w) <C 1 by assumption. 

The term 1 — 2Z 2 in equation ()35j) is due to the fact that we express the equation of 
motion of S with respect to the Fermi time T rather than r, the proper time of S. For 
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FIG. 2: Plot of f versus T/(GM) based on the integration of equation 1)351) for r$ = 5GM, p = 2 
in the top panel and p = —2 in the bottom panel, a = GM and Tq = \/2, 2, 3, 4 and 5. 

\Z\ <C 1, equation (|Ho]) reduces to a Jacobi-type equation for the relative motion of S with 
respect to O 

To characterize the relative motion invariantly, let us imagine a set of static observers in 
the Fermi system situated along the Z axis and let T be the Lorentz factor of S as measured 
by these static Fermi observers that are in general accelerated. Then f = y/— goo r, where 
—goo = 1 + kZ 2 in our approximation scheme; therefore, T < T, since k < 0. 

Equation contains the "electric" and "magnetic" curvatures k and q, respectively, 
that can be expressed as functions of the Fermi coordinate time T. That is, the equation 
for r in ()32|) can be integrated such that with s i— > T, we find r(T) and hence fc(T) and 
q(T) that must then be substituted in (|35|) . This system has been numerically integrated 
with initial conditions that at T = 0, Z = and Z = Vq > 0, and the results are presented 
in figures 1 and 2, where T is plotted versus T/(GM) for ultrarelativistic relative motion 
along the rotation axis of a maximal (a = GM) Kerr source. We note that with our initial 
conditions both T and T at T = are given by Tq = (1 — Vq 2 ) -1 / 2 , which is the Lorentz 
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factor corresponding to Vq < 1. For motion along the Z axis, T < T; however, it turns 
out that a similar plot for F would be indistinguishable from figure ^ We choose 7 = 1 in 
figures 1 and 2, since the results turn out to be independent of the choice of 7 so long as O 
is slowly outgoing. For highly ultrarelativistic motion as in figure ^ the spin of the clump 
has a negligible influence on its deceleration. In fact, with the same conditions as in figure Q 
but for p = and 10, the resulting figures turn out to be indistinguishable from figure d 
Let us note that r — > 00 is not allowed here, since the back reaction of such a particle 
on the gravitational field of the source cannot be neglected and hence our test particle 
approximation scheme would break down. The case where the spin of S is antiparallel to 
the Z axis can be treated by formally letting s l— > ~ s o- The corresponding results are 
presented in the bottom panel of figure El The graphs in the top panel are slightly above 
the corresponding ones in the bottom panel, though the trends are essentially the same. 

Let us suppose that the clump is homogeneous and cylindrical with radius R, and the 
axis of the cylinder coincides with the Z axis. Thus so = ^mR 2 u, where to is the frequency 
of rotation of the clump such that Ruo < 1. Moreover, (GM/r)p = so/(mr) <C 1 implies that 
p <C r/(GM). In figures H and |2] we have r > ro = 5GM; therefore, we have chosen p = 2 
to represent in effect the maximum spin of the clump. This choice is consistent with various 
scenarios that have been proposed for jet creation near the poles of a rapidly rotating Kerr 
black hole (see [lf|, chapters 9 and 11, and the references cited therein). 



VI. DISCUSSION 

If the initial Lorentz factor r of S is close to the critical value T c = v2, then the spin- 
curvature force is rather small but may not be negligible very close to the central black 
hole. However, for r ^> y/2, the contribution of the spin-curvature force turns out to be 
essentially negligible and the spinning particle S decelerates toward the critical speed. This 
confirms our previous work on the speed of jets in microquasars ~v\ . 

In this paper, we have shown that the current speculation that the spin-curvature cou- 
pling could be a prominent driver of astrophysical jets from black holes is overly optimistic. 



We show quite generally that the standard tidal effects 



3 wi 



will dominate the spin-curvature 



11 



effects for plasma clumps in astrophysical jets. 
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